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Abstract 

We compute the dynamical prefactor in the nucleation rate of bubbles or droplets 
in first order phase transitions for the case where both viscous damping and ther- 
mal dissipation are significant. This result, which generalizes previous work on 
nucleation, may be applied to study the growth of bubbles or droplets in condensed 
matter systems as well as in heavy ion collisions and in the expansion of the early 
universe. 
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Near the critical temperature for a first order phase transition, matter in a 
certain phase is metastable and small fluctuations in the state variables may acti- 
vate the nucleation of bubbles (a term we use below to represent either bubbles or 
droplets) of a more stable phase of matter. If these bubbles have a radius larger than 
some critical value, they begin to grow exponentially The nucleation and growth 
of these critical bubbles has historically been of great interest in the physics of 
liquid-gas phase transitions and in condensed matter physics Q. More recently, 
they have been studied in the context of first order phase transitions in the early 
universe |2|, |^, ^ and in high energy heavy ion collisions |5| . 

A general kinetic theory of homogeneous nucleation was developed by Langer Q . 
In the neighbourhood of a first order phase transition, when the critical radii of the 
bubbles exceed the correlation length, a reduced description of nucleation in terms 
of a coarse-grained free energy is appropriate. Langer and Turski used such a 
phenomenological approach to show that the nucleation rate of bubbles could be 
written as a product of three terms 

/= ^J7oexp(-AF/r). (1) 

The dynamical prefactor k determines the exponential rate of growth of bubbles of 
a critical radius, beyond which size the bubbles are stable. For the bubbles to grow 
beyond the critical radius, latent heat must be carried away from the surface of 
the bubble. This is achieved through thermal dissipation and/or viscous damping. 
Kawasaki and Turski and Langer , neglecting viscous damping, show that k for 
the condensation of a supersaturated vapour is linearly proportional to the thermal 
conductivity of the vapour. The statistical prefactor 0,q is a measure of the phase 
space volume of the saddle point region of the free energy functional and AF is the 
change in the free energy required to activate the formation of a critical bubble. 

Recently, the theory of Langer and Turski has been used by Csernai and Ka- 
pusta to study nucleation in relativistic first-order phase transitions Q . The baryon 
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density in the systems studied is negligible. In the absence of a net conserved charge, 
the thermal conductivity vanishes and the expression of Langer and Turski for the 
dynamical prefactor is no longer applicable. It was shown in Ref. that, for the 
systems studied, the latent heat could be transported from the growing bubble by 
viscous damping instead of thermal dissipation; the new expression for the dynam- 
ical prefactor depends linearly on the shear viscosity of the surrounding medium. 
The resulting expression for the pre-exponential factor differs significantly from ear- 
lier estimates where, on dimensional grounds, the prefactor was taken to be or 
T^. Here T is the temperature and Tc is the critical temperature of the first order 
phase transition. 

In this work, we derive a general expression for the dynamical prefactor in the 
nucleation rate of critical sized bubbles in first order phase transitions. This formula 
may be used to study nucleation in liquids and gases and in condensed matter 
systems where both the viscosity and thermal conductivity are significant. It may 
also be used to estimate the probability of formation of superheated quark-gluon 
droplets or supercooled hadronic bubbles in the baryon-rich matter produced in high 
energy heavy ion collisions Further applications also include the expansion of 
the early universe and the formation of neutron stars. 

Consider two phases of matter, phase A and phase B, where phase ^ is a 
metastable state which decays into the more stable phase B. Phase A, for instance, 
may be a supersaturated vapour which when supercooled, nucleates bubbles of a 
liquid phase B. When the radius i?* of the critical bubble of phase B is greater than 
the correlation length ^ in phase A, the behaviour of the system can be described in 
terms of a coarse-grained energy functional F. This functional depends on the fiuid 
density n(r), the local fiuid velocity U{r) and the temperature T{r). Following 
Ref. [^, we make the ansatz that F = Fk + Fi where Fk is the kinetic energy. 
The interaction term Fi is the sum of the Helmholtz free energy and the van der 
Waals-Cahn-Hilliard gradient energy The formalism for the coarse-grained 
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free energy functional is also valid for relativistic systems where the pressure is 
comparable to the energy density and the fluid velocity U{r) in the local rest frame 
is small compared to the speed of light. 

In addition to the homogeneous phases A and B, the free energy is also station- 
ary for a configuration whose solution is a generalization of the van der Waals soliton, 
which has the hyperbolic tangent -like density profile. To determine the expansion 
of the bubble about this stationary configuration, we linearize the hydrodynamic 
equations around the stationary configuration: n{f) = n[r) + v[r)^ U{r) = + U{f) 
and T(r) = Tq + 9{f), where the quantities z^, U and 9 correspond to small devi- 
ations in the density, velocity and temperature, respectively, from their stationary 
values. They approach constant values away from the interface. We derive below a 
general expression for k. which does not depend on any specific parametrization of 
the free energy. 

In Ref. @ relations were derived between the velocity potential and the density 
and temperature functions on either side of a diffuse interface. These relations, the 



Kotchine conditions [12|, are generalizations of the well-known Rankine-Hugoniot 
discontinuity conditions for shocks. For instance, in the former case, the velocity 
of matter diffusing through the interface is a function of position and falls off away 
from the interface. In the later case, the velocity of the matter is a constant. It 
was shown in Ref. Q that these Kotchine conditions give the correct dispersion 
spectrum for capillary waves. 

We now use the Kotchine conditions for a spherically growing bubble to derive 
an expression for the dynamical prefactor k. Our derivation is similar to that 
of Turski and Langer but differs from theirs in some key aspects. The Kotchine 
conditions for a spherical bubble are 

[nUR] = [nQ, (2) 



M=0, (4) 

In the above, the brackets denote the difference in the bracketed quantity across the 
interface. For instance, [n] = ni, — Ua = An, where the subscripts denote phase B 
and phase A, respectively. Also, Ur is the velocity of matter through the interface, 
dR/dt is the velocity of the bubble wall, P is the pressure and ^ the chemical 
potential. The latent heat per particle is given by /, A is the thermal conductivity 
and rj and C, are the shear and bulk viscosities respectively. 

The first Kotchine condition, Eq. (^), is the matter continuity relation across 
the interface. The second Kotchine condition is the well-known Laplace formula for 
the surface tension. The third denotes the continuity of chemical potentials between 
the two phases at Tc- The final Kotchine condition equates the latent heat produced 
per unit area per unit time at the interface to the energy dissipated per unit area 
per unit time. 

Combining Eq. (|2|) and Eq. (|5|), the total energy flux transported outwards is 
given by 

dR dT (A \ dUR 

Here Aw is the difference in the enthalpy densities of the two phases. From the 
continuity relation, dti' = —V • {fiUfi), one may show (see the discussion preceding 
Eq. (77) in Ref. Q) on very general grounds that the radial dependence of the veloc- 
ity at the interface Ur oc Hence, dU{r)/dr\r=R = — 2C/^(r)/i?. Substituting 
this relation in the above equation, we obtain 

^ dR dT ^ /4 A Ul 



dt dr V3 J R 

We wish to obtain a similarly simple expression for the gradient in the temper- 
ature dT/dr. For the systems we consider, we may assume that the temperature 



varies slowly across the bubble wall. If we represent the temperatures in the two 
phases by Ta and Tf,, we can define an average temperature T and the variation 
hy Ta = T + and = T — 6q. In the quasi-stationary approximation ~ 0. 
The solution to the Laplacian is then ||l^ 

9b = Oo yr<R, 



yr>R, (8) 



where Oq is a constant. Hence, 



-I --^ (9) 
dr K=R R ■ 

To determine we use the continuity of chemical potentials /x^ = across the 
interface of the critical bubble (the third Kotchine condition). Then, using the 
first law of thermodynamics and assuming a large latent heat (strong first order 
transition), we arrive at the relation 

^-^--f. (10) 

nb Ua T 
Now from the second Kotchine condition, Eq. (^), 

n = Pa + |^. (11) 

Substituting this equation in Eq. (|To|), we obtain, after a little algebra 

Pa = ^f- + ^V (12) 



An\R T , 

The temperature difference between the two phases, 29q, is due to the dissipation 
of latent heat. For the critical bubble, Bq = 0, which implies that 

Replacing Pa in Eq. (^) with the above expression, we obtain finally for 6o the 
relation 

2aT ( 1 1 



AwR \R^ R 



(14) 



Substituting this result for 6*0 in Eq. (^), we have 

dT, 2aT ( 1 1 



(15) 



dr 'r AwR V-R* 
We have one further unknown-the velocity Uii{r) of matter diffusing through 
the surface of the growing bubble. If there exists a net momentum flux through the 
interface, then from Laplace's formula, 

Au;t/i = 2a(^i--l) . (16) 

We have omitted the shear term in the above equation since it represents a higher 
order contribution to the linearized hydrodynamic equations. 

Combining our results in Eq. (|l|) and Eq. (||) with Eq. (0), we obtain the 
following expression 



dR 2g ( \ 1 



dt {AwYR \R^ R 



AT + 2 ( ^r/ + C 



(17) 



\iR — R^ oc exp(Kt), we obtain finally our general result for the dynamical prefactor 
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{Aw^Rl 



(18) 



In the limit of zero baryon number, A ^ 0, and we obtain the result of Csernai and 
Kapusta. If the matter is baryon-rich but viscous damping is negligible, 77, (" — > 0, 
we obtain the result of Kawasaki, and Turski and Langer. Since the results for 
the dynamical prefactor were known in the two limits, and since only terms linear 
in the transport coefficients are retained, the above result might easily have been 
intuited. However, many of the ingredients in our calculation are generally valid, 
and hence the dynamical prefactor could also be computed in a similar manner for 
more complicated systems. 

We should point out that there are several assumptions that have been made 
in our derivation of the dynamical prefactor. Our result is strictly valid when 
non-linear effects can be ignored and the linearized hydrodynamic equations are 



applicable. Further, for the coarse-graining description to hold, the radii of the 
bubbles must be larger than the correlation length. We have also assumed that 
heating due to dissipation is slow, causing the temperature to vary slowly across 
the bubble wall. Finally, we have assumed in our derivation that the phase transition 
is strongly first order-releasing considerable latent heat. 

To summarize, we have derived above an expression for the dynamical prefactor 
which governs the growth of critically sized bubbles nucleated in first order phase 
transitions. Our results are applicable to the wide range of phenomena where both 
viscous damping and thermal dissipation effects are important. In a following pa- 
per we will discuss one such application-the nucleation of quark-gluon droplets 
in baryon-rich hadronic matter created in high energy heavy ion collisions. 
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